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ABSTRACT 


This  report  is  the  last  of  a  series  of  three  reports  on  a  com¬ 
prehensive  study  of  CTD  instrument  lowering  mechanics.  The  first 
report,  WHOI  79-81,  "A  Study  of  CTD  Cables  and  Lowering  Systems”, 
examines  the  causes  and  modes  of  lowering  cable  failures,  both  mechanical 
and  electrical,  and  makes  recommendations  to  improve  existing  instrument 
packages  and  lowering  procedures.  The  second  report,  WHOI  81-76, 
"Hydrodynamics  of  CTD  Instrument  Packages",  is  a  theoretical  study  of 
instrument  package  stability  when  cable  lowered  or  free  falling.  The 
model  is  used  to  predict  the  hydrodynamic  response  of  CTD  packages  in 
their  present  or  improved  configuration.  This  report,  WHOI  83-21,  is 
more  factual.  It  describes  the  tests  performed  on  scale  models  and 
actual  CTD  packages  to  actually  observe  and/or  measure  their  hydrodynamic 
behavior.  Analytical  results  and  experimental  data  obtained  in  this 
study  are  used  to  draw  recommendations  for  CTD  package  improvement  and 
future  lowering  procedures. 
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1.0  BACKGROUND 

Frequent  failures  of  cables  used  for  lowering  instrument 
packages  from  vessels  at  sea  (Figure  1)  have  prompted  a  study  to  identify 
the  causes  of  these  failures  and  to  recommend  improvements  in  materials, 
handling  equipment  and  procedures.  While  this  study  was  specifically 
“"cdhcetned  with  Conductivity,  Temperature  and  Depth  (CTD)  packages  it  has 
relevance  to  many  cable  lowered  packages.  A  typical  CTD  instrument 
package  is  shown  in  Figure  2. 

The  first  part  of  the  study  (Berteaux  et  al^  1979)  examined  the 
causes  and  modes  of  cable  failure  both  mechanical  and  electrical.  Fail¬ 
ures  as  a  consequence  of  high  tensile  loads,  fatigue,  corrosion,  mishand¬ 
ling,  shock  loads,  cross  winding  and  jumped  sheaves  were  identified.  Methods 
to  improve  the  reliability  of  CTD  cables  and  lowering  systems  were  proposed. 
This  included  quality  control  and  acceptance  testing  procedures  for  new 
cables  and  instructions  for  installing  these  cables  aboard  ship.  Proper 
maintenance  procedures  for  both  cables  and  winches  were  also  reviewed. 
Operational  limits  for  both  depth  of  cast  and  payout  speed  were  calculated 
for  different  sea  states.  Alternate  CTD  lowering  systems  using  special  steel 
armored  cables,  synthetic  fiber  and  fiber  optics  cables  were  evaluated. 

The  results  of  this  investigation  were  condensed  in  a  first  report 
(Reference  1).  The  report  concludes  with  recommendations  including  a  re¬ 
design  of  the  CTD  instrument  package  to  improve  its  hydrodynamic  performance 
characteristics. 

A  second  report  (Cook  1981,  Reference  2)  studied  the  hydrodynamics  of 
CTD  instrument  packages.  In  this  report  an  analysis  of  the  standard  Woods 
Hole  CTD  package  was  made.  Its  static  stability  was  investigated  with  the 
help  of  a  computer  model.  The  model  was  used  to  perform  sensitivity 
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analyses  considering  modifications  to  the  standard  package  in  symmetry, 
weight  and  drag.  From  thia,  an  optimal  configuration  was  selected. 

Criteria  for  the  formulation  of  a  half  scale  model  were  developed  to 
satisfy  geometric,  kinematic  and  dynamic  similarity.  The  report  concluded 
with  recommendations  to  improve  the  hydrodynamic  behavior  of  the  package 
through  an  increase  in  its  terminal  velocity  by  changes  in  symmetry,  drag 
reduction  and  the  addition  of  weight. 

A  paper  reviewing  these  studies  and  presenting  the  results  of  tank 
model  tests  was  published  in  Deep  Sea  Research  (Berteaux  and  Walden,  1983, 
Reference  3).  This  paper  also  presents  the  results  obtained  at  sea  which 
confirmed  our  suspicions  that  conditions  of  zero  tension  in  the  wire  did 
indeed  occur  when  the  combination  of  ship  motions  and  payout  speed  were  too 
great.  Severe  wire  damage  could  occur  from  the  subsequent  snap  loading. 

The  report  which  follows  concludes  our  CTD  study.  It  presents  in 
detail  the  instruments  and  procedures  used  to  observe  and  document  the 
hydrodynamic  behavior  of  standard  and  improved  CTD  packages.  This  report 
also  presents  a  number  of  recommendations  to  improve  the  efficiency  and 
the  reliability  of  instrxraient  lowerings. 

The  first  section  describes  how  the  terminal  velocity  of  free  falling 
CTD  packages  was  actually  measured.  The  effects  that  drag  reduction  and 
increased  weight  have  on  package  terminal  velocity  are  discussed  and 
numerical  results  presented. 

A  half  scale  model  of  the  CTD  package  was  built  and  its  flight 
pattern  observed  as  it  was  dropped  in  the  large  water  tank  of  the  Naval 
Surface  Weapons  Center,  The  second  section  describes  the  criteria  of 
similarity  used  for  the  fabrication  of  the  model;  it  outlines  the  test 
procedures,  and  reviews  the  results  obtained. 

The  next  section  covers  the  tests  performed  at  sea  during  actual  CTD 
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f  lowerings.  It  first  describes  the  tests  done  with  a  standard  CTD  package. 

This  package  was  instrumented  to  prnvide  simultaneous  measurements  of 
■  package  depth,  inclination  arid  cable  tension  at  the  point  of  package 
r attachment.  Important  results  obtained  with  these  measurements  are 
;  reviewed.  Comments  on  the  performance  at  sea  of  an  improved  compact  CTD 
package  are  next  presented. 

In  the  final  section  the  report  formulates  a  number  of  recommenda¬ 
tions  to  improve  the  efficiency  and  the  reliability  of  instrument  package 
lowering  operations.  Prudence  and  common  sense  will  dictate  that  limits 
be  set  on  payout  speeds  and  lengths  payed  out  as  a  function  of  sea  state. 
Theoretical  considerations  and  measurements  made  during  the  study  permit 
quantifying  these  limits.  Ways  of  specifically  improving  CTD  packages 
are  discussed.  A  brief  review  of  ship  motion  compensation  and  of  its 
benefits  for  cable  lowering  applications  concludes  this  report. 

2.0  ACTUAL  MEASUREMENTS  OF  CTD  PACKAGE  TERMINAL  VELOCITY 

2.1  Instrimentation 

An  instrument  was  constructed  to  measure  actual  velocities 
of  CTD  packages  in  shallow  water  dock  tests  (Figure  3) .  The  instrument 
consisted  of  a  lightweight  bicycle  wheel  rim  geared  through  a  friction 
clutch  to  a  fifteen  turn  linear  potentiometer.  Small  Kevlar  line 
(.045  inch  OD)  was  wound  on  the  rim  and  attached  to  the  package  to  be 
dropped.  The  potentiometer  was  connected  to  a  DC  power  supply  and  to  a 
strip  chart  recorder  which  displayed  voltage  (distance  in  this  case)  versus 
time.  By  differentiating,  or  taking  the  slope  of  the  distance  versus 
time  curve,  velocity  as  a  function  of  time  was  obtained.  Terminal  velocity 
can  be  determined  from  the  constant  slope  section  of  the  distance  versus 
time  plot. 
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FIGURE  3.  Instrument  to  measure  terminal  velocity 


2.2  Procedure 


The  operating  technique  for  using  this  apparatus  was  to 
first  zero  the  recorder  with  the  CTD  just  below  the  surface.  Then  ad¬ 
justments  were  made  in  the  power  supply  output  until  full  scale  reading 
was  obtained  with  the  package  near  the  bottom.  At  the  dock  the  water 
depth  is  approximately  64  feet.  The  recorder  was  turned  on  just  prior  to 
quick  releasing  of  the  package.  As  the  CTD  fell,  the  instrument  measured 
and  recorded  the  distance  traveled  versus  time  until  a  restraining  line 
brought  the  package  to  rest  a  few  feet  off  the  bottom.  At  this  time  the 
recorder  was  shut  down. 

2.3  Results 

Early  velocity  measurements  were  made  of  the  standard  WHOI 
CTD  package  depicted  in  Figure  4.  This  package  consisted  of  a  CTD  instru¬ 
ment  with  a  24  inch  pressure  case,  a  (24)  1.2  liter  bottle  rosette  sampler, 
and  a  pinger  all  enclosed  in  a  3/4  inch  galvanized  steel  pipe  frame.  In 
water  the  90  pound  CTD  instrument  makes  up  35%  of  the  total  package  immersed 
weight  (254  pounds),  the  70  pound  rosette  287o,  the  40  pound  pinger  167o 
and  the  54  pound  frame  the  remaining  21%.  Lead  pigs  clamped  to  the  lower 
part  of  the  CTD  instrument  case,  were  added  and  changes  iji  terminal  velocity 
of  the  standard  package  due  to  increased  weight  were  measured  and  recorded. 

The  standard  package  was  then  reconfigured  and  ballasted  to  obtain  a 
higher  terminal  velocity  with  less  weight.  The  external  steel  frame  was 
completely  removed  except  for  the  solid  rod,  stainless  steel  frame  at  the 
bottom  protecting  the  CTD  sensors.  This  substantially  decreased  the  package 
hydrodynamic  drag.  The  rosette  housing  was  attached  closely  to  the  upper 
CTD  end  cap  by  stand  offs  further  reducing  drag.  Weight  was  added  to  this 
compact  package  by  clamping  a  smooth  lead  sleeve  to  the  lower  end  of  the 
CTD  case,  as  shown  in  Figure  5, 
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FIGURE  5 
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Records  of  distance  (depth)  versus  time  for  a  standard  and  a 
compact  CTD  package  are  shovm  in  Figure  6.  The  full  scale  depth  in 
this  plot  is  63.75  feet  with  a  chart  speed  of  one  inch  per  second.  The 
immersed  weight  of  the  standard  package  was  255  pounds  with  an  air  weight 
of  420  pounds.  It  reached  a  terminal  velocity  of  7.59  ft/sec  in  about 
two  seconds  after  drop.  A  similar  trace  is  shown  for  the  compact  version 
having  an  immersed  weight  of  292  pounds  and  an  air  weight  of  434  pounds. 

In  this  case  a  terminal  velocity  of  10.8  ft/sec  was  reached  in  just  over 
one  second. 

Table  1  lists  the  terminal  velocity  for  five  configurations,  three 
for  the  standard  version  and  two  for  the  compact  configuration.  Different 
ballast  weights  are  shown  for  each  version.  The  results  clearly  show 
the  effect  of  increasing  weight  and  decreasing  area  on  the  terminal 
velocity.  By  adding  lead  weight  in  104  pound  wet  weight  increments  to  the 
standard  package  the  7.59  ft/sec  velocity  was  increased  to  9.04  ft/ sec 
and  finally  to  10.31  ft/sec  with  207  pounds  of  lead  added.  This  increased 
the  terminal  velocity  by  377®,  but  also  increased  the  weight  to  648  lbs 
in  air. 

The  compact  package,  CTD  and  rosette  with  no  external  frame,  initially 
had  a  velocity  of  8.69  ft/sec  with  a  wet  weight  of  only  125  pounds.  By 
adding  a  lead  sleeve  of  167  lbs  (wet  weight) ,  the  velocity  of  this  compact 
package  was  increased  to  10.8  ft/sec.  The  immersed  weight  of  this  last 
version (292  lbs)  is  almost  equal  to  the  immersed  weight  of  the  standard 
package  (255  lbs).  Yet,  because  of  its  compact  shape,  its  terminal 
velocity  is  43%  larger.  Its  weight  in  air  is  only  434  lbs,  which  is  about 
the  same  as  the  420  lbs  air  weight  of  the  standard  package. 


FIGURE  6.  Curves  of  depth  versus  time  of  two  free  falling  CTD  packages 
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Terminal  Velocities 
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255 

420 

7.59 

- 
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+(4)  lead  weights 

359 
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9.04 

20% 

Standard 

+(8)  lead  weights 
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648 

10.31 

37% 

Compact 

125 

240 

8.69 

15% 

Compact 

with  lead  sleeve 

292 

434 

10.8 

43% 

Reducing  the  package  frontal  area  also  reduces  the  amount  of  entrained 
water  thus  improving  the  package  ability  to  accelerate  and  to  better 
accomodate  ship  roll  and  heave.  The  large  righting  moment  resulting  from 
the  drag  forces  on  the  rosette  bottles  located  high  above  the  center  of 
gravity  combined  with  the  gravity  pull  of  the  lead  sleeve  low  on  the  CTD 
produces  very  stable  flight  characteristics. 

3.0  SCA.LE  MODEL  TESTS 

A  series  of  CTD  scale  model  tests  were  conducted  in  July  1981 
at  the  Naval  Surface  Weapons  Center  (NSWC) ,  Silver  Spring,  Maryland.  The 
objective  of  these  tests  was  to  observe  and  record  the  "flight”  pattern  of 
scaled  down  versions  of  CTD  packages  as  they  were  cable  lowered  and/or 
allowed  to  free  fall  to  the  bottom  of  the  NSCW  vertical  tank.  This  tank  is 


X5  meters  in  diameter  and  30  meters  deep. 

3.1  Criteria  of  Flow  Similarity 

The  hydrodynamic  behavior  of  a  full  scale  prototype  can 
be  reasonably  assessed  from  tests  performed  with  a  scale  model  provided 
certain  conditions  of  similarity  between  prototype  and  model  are  repro¬ 
duced.  In  short,  one  must  ensure  that: 

1.  The  test  procedure  introduces  the  same  effects  on  the  model. 

2.  These  effects  cause  the  same  results. 

Criteria  of  similarity  commonly  used  include: 

o  Geometric  similarity.  The  model  must  be  a  good  physical  image 
of  the  prototype.  It  must  reproduce  its  shape  and  also  its 
main  physical  characteristics.  For  example  identical  location 
of  center  of  buoyancy,  center  of  gravity,  similar  mass  moments 
of  inertia,  etc... 

o  Kinematic  similarity  which  ensures  that  the  streamline  pattern 
remains  the  same  for  model  and  prototype. 

o  Dynamic  similarity,  which  ensures  that  model  and  prototypes 

are  subjected  to  the  same  loads.  This  is  accomplished  by  main¬ 
taining  the  ratios  of  homologous  forces  identical  for  the 
model  and  the  prototype. 

Given  the  size  of  the  actual  CTD  package,  building  a  half  scale 
model  seemed  reasonable  and  practical.  Its  weight  was  determined  using 
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I  According  to  the  criterion  of  dynamic  similarity,  the  ratios  of  the 
avity  and  drag  forces  applied  to  the  model  and  package  had  to  be 
e  same.  Using  the  familiar  expression  for  the  drag  force 


Ft>-  “i:  ^  Cp  Av  ^ 

where  ^  is  the  fluid  density 

is  the  drag  coefficient 
A  is  the  cross  section  area  given  by 
^  being  some  geometric  constant  (i.e. 
and  D  a  characteristic  dimension 
V  the  speed  of  fall 
the  ratios  equality  can  be  written  as: 


<\  D  ^ 


CopSpOpWp 


(3.1) 


where  W  is  the  immersed  weight,  and  the  subscripts  M  and  P  refer  to 
model  and  package  respectively. 

In  flow  regimes  dominated  by  inertia  and  viscous  forces,  dynamic  similarity 
will  also  dictate  that  the  Reynold *s  number  be  the  same  for  both  model  and 
prototype.  Both  being  immersed  in  water,  an  equality  of  Reynold's  numbers 


yields: 


V  o  =  V  D 

M%\ 


(3.2) 


For  geometric  similarity  to  prevail 

Saa  ~  S  p 


(3.3) 


Model  and  package  having  the  same  Reynold's  numbers  and  the  same  shape, 
both  have  the  same  drag  coefficient,  thus 

==  Ce»c. 


(3.4) 


Introducing  equalities  3.2,  3.3,  and  3.4  in  3.1  yields 


=  Wp 

Thus  the  model  immersed  weight  should  equal  the  inmersed  weight  of  the 
prototype.  Introducing  the  scale  ratio  =  2  in  3.2  yields 

V  ^  2  V 

A  half  scale  model  which  falls  twice  as  fast  and  weighs  in  water  as  much 
as  the  actual  package  is  not  a  very  practical  one.  A  lighter  model 
traveling  at  moderate  speeds  would  be  easier  to  handle  and  to  visually 
observe.  In  fact  such  a  model  may  not  be  as  distorted  as  it  may  seem. 
Certainly  the  flow  pattern  around  the  blunt  surfaces  of  the  package  compon¬ 
ents  is  turbulent  most  of  the  time.  This  implies  that  the  drag  coefficient 
remains  constant  over  a  large  range  of  Reynold's  numbers.  If  the  speed 
of  the  model  is  kept  large  enough  to  maintain  turbulent  flow  then  its  own 
drag  coefficient  will  equal  the  drag  coefficient  of  the  prototype  over  the 
same  range  of  Reynold's  numbers.  In  this  case  the  criterion  of  dynamic 
similarity  3.1  reduces  to 


w. 


-  ^ 


from  which 

~  ^  \  Vp  /  (3.5) 

The  package  terminal  velocity  is  of  the  order  of  10  ft/sec.  Selecting  a 


speed  of  model  fall  equal  to  the  speed  of  prototype  fall  -  which  should  be 
sufficient  to  maintain  turbulent  flow  conditions  -  yields 

w  = 

M  ^ 


To  summarize,  the  model  retained  would  be  half  scale  and  would  have  a 
submerged  weight  equal  to  one  quarter  of  the  full  scale  value.  At  terminal 
velocity,  the  model  will  have  a  speed  given  by 


V  - 

^  yCoM  yCo^  AV/4* 


which  equals  the  terminal  velocity  of  the  package  as  long  as  =  Cj^p 

3.2  Design  and  Construction  of  Scale  Model 

The  half  scale  model  was  made  of  removable  cylinders 
attached  to  a  tubular  frame  (Figure  7) .  These  cylinders  had  half  the 
diameter  and  half  the  length  of  the  components  represented.  They  were 
placed  in  the  frame  at  a  location  similar  to  the  location  occupied  by  the 
component  in  the  actual  package.  To  maintain  correct  submerged  weight 
distribution  as  well  as  correct  weight-to-drag  ratios,  the  material  density 
of  each  cylinder  was  calculated  according  to  the  following  scheme: 

1.  Measure  the  immersed  weight  of  the  actual  component  "i",  Wp  f  t ) 

2.  The  immersed  weight  of  the  cylinder  representing  the  component 
’’i"  is  then 

'A/«vf.i)  *•  WpfO/oil. 

3.  Compute  the  volume  VC)^(i)  of  the  cylinder  ”i” 

4.  The  desired  specific  weight  of  the  cylinder  "i"  is  then  given  by: 

THe.  oeusi-rr  of  seA 


Results  obtained  following  this  approach  are  summarized  in  Table  2, 


TABLE  2 


Specifications  for  \  Scale  CTD  Model 


Component 

Modeled 

Cylinder  Dimensions 
(inches) 

Diam  x  Length 

Weight 

(lbs) 

Immersed 

Weight 

(lbs) 

Specific 
Weight 
(Ibs/cu/ ft) 

Specific 

Gravity 

CTD 

3.5  X  12.5 

29.26 

24.75 

415.63 

6.6 

Rosette 

11.5  X  15.0 

65.20 

17.50 

83.41 

1.34 

Frame 

\  0  diam. tubing 
Ring  &  Verticals 
to  scale 

16.90 

13.50 

310.0 

4.97 

Pinger 

2.5  X  13 

12.24 

10.0 

349.71 

5.60 

Nephelometer 

3.0  X  16 

22.87 

18.68 

349.71 

5.60 

The  immersed  weight  of  an  actual  CTD  package  made  up  of  a  CTD  instrument,  a 
rosette,  a  WHOI  frame  and  one  pinger  was  measured  to  be  255  lbs.  The 
immersed  weight  of  the  model,  using  the  data  shown  in  the  table,  is  found 
to  be  65.75  lbs.  or  3.9  as  small  as  the  prototype. 

Plates  of  steel,  aluminum,  and  polyethylene  were  used  to  fabricate  the 
model.  Circular  slabs  of  these  materials  were  stacked  in  combinations 
yielding  the  proper  dimensions,  immersed  weight,  and  location  of  center  of 
gravity  for  each  cylinder. 

3.3  Model  Test  Procedure 

Models  of  instrument  packages  were  lowered  with  the  help 
of  an  adjustable  speed  winch.  Maximum  winch  speed  exceeded  the  terminal 
velocity  of  all  model  configurations  tested  thus  allowing  a  number  of 
lowerings  to  be  essentially  free  fall. 


The  lowering  cable  was  made  somewhat  shorter  than  the  tank  depth. 

The  models  were  attached  to  the  lowering  cable  by  a  length  of  compliant 
nylon  rope  which  could  absorb  the  kinetic  energy  of  the  lowered  or  free 
falling  specimen  at  the  end  of  payout.  Keeping  the  winch  running  in  the  same 
(payout)  direction  simply  brought  the  specimen  back  to  the  surface. 

The  cable  was  passed  through  an  overhead  sheave  located  above  the 
tank  center.  A  load  cell  attached  to  the  sheave  was  used  to  monitor  and 
record  tension  during  the  lowerings. 

All  lowerings  were  made  a  constant  speed.  The  lowering  speed  was 
inferred  from  rpm  measurements  made  with  a  digital  tachometer.  Typical 
lowering  speeds  were  4.3,  6.7,  and  10  ft/sec,  (78.6,  122.5  and 
182.9  meters/min. ) . 

Flight  attitude  and  deviation  from  vertical  plumb  line  were  observed, 
noted,  and  in  certain  cases  filmed  with  the  help  of  horizontally  and 
vertically  mounted  underwater  cameras.  To  help  visualize  the  flight  path 
as  seen  from  the  top,  two  yellow  polypropylene  lines  crossing  each  other 
at  right  angles  were  tied  to  the  tank  bottom  grid. 

Various  combinations  of  instrument  models,  stabilizing  fins  and  nose 
cone  were  used.  The  different  configurations  tested  included: 
o  CTD,  Rosette,  used  as  "base  line”, 
o  CTD,  Rosette,  Finger  strapped  to  the  frame, 
o  CTD,  Rosette,  Finger,  Nephelometer  positioned  horizontally 
across  the  bottom  of  the  frame,  on  center  and  off  center, 
o  CTD,  Rosette,  Nephelometer  both  on  and  off  center, 
o  CTD,  Rosette,  3  small  vertical  fins, 

o  Same  as  above  with  3  large  vertical  fins. 

o  Same  as  above  with  addition  of  pinger  strapped  on  the  side, 

o  CTD,  Rosette,  Finger,  Nose  Cone,  3  large  fins. 


Test  Results  and  Conclusions 


The  following  comments  summarize  the  test  results: 

Terminal  velocities  observed  agree  with  terminal  velocities 
predicted. 

At  small  lowering  speeds  and  therefore  high  cable  tension  values 
all  configurations  behave  rather  well.  Conversely  deviations 
from  well  behaved  plumb  line  flight  could  be  observed  only  at 
high  lowering  speeds  and/or  at  terminal  velocities  (free  fall). 
At  speeds  of  6.0  f t/sec  up  to  terminal  velocity  (7.6  ft/sec) 
the  baseline  configuration  (frame,  rosette  and  CTD  instrument 
casing)  is  statically  stable  but  dynamically  unstable.  It 
falls  in  a  straight  line  but  oscillates  back  and  forth  around 
its  c.g.  .  The  pitch  amplitude  appears  to  be  approximately 
one  CTD  diameter  and  the  frequency  one  hertz. 

Strapping  a  pinger  outside  the  frame  has  the  following  effects; 

It  makes  the  package  tilt  about  ten  degrees. 

The  added  weight  causes  an  increase  in  terminal  velocity. 

The  added  weight  lowers  the  package  center  of  gravity 
and  increases  the  distance  between  package  center  of  gravity 
and  center  of  buoyancy.  As  a  result  the  wobbling  of  the 
baseline  configuration  is  sensibly  reduced  but  still  present. 
Strapping  a  nephelometer  across  the  frame  bottom  ring  produces 
the  following  effects: 

The  added  weight  causes  an  increase  in  terminal  velocity. 

The  location  of  the  added  weight  on  the  frame  bottom  consid¬ 
erably  increases  the  distance  between  package  center  of 
buoyancy  and  center  of  gravity.  As  a  result  the  oscillations 


observed  in  the  baseline  are  totally  suppressed. 

Horizontally  mounted  cylinders  seem  to  induce  spin.  They 
make  the  package  turn  around  its  vertical  axis  as  it  travels 
through  water,  more  turns  being  induced  at  higher  speed 
(typically  two  turns  in  100  feet  drop  at  terminal  velocity). 
The  addition  of  a  nose  cone  did  produce  some  unexpected  (and 
spectacular)  results. 

It  did  not  cause  an  increase  in  terminal  velocity. 

It  produced  static  instability  resulting  in  a  large  radius 
helical  flight  path. 

The  addition  of  three  equally  spaced  vertical  stabilizing  fins 
resulting  in 

Total  suppression  of  baseline  oscillations. 

Absolute  pltimb  line  vertical  flight  pattern  for  concentric 
(axl symmetrical)  configurations. 

With  a  pinger  attached  on  the  side  of  the  frame  the  package 
will  fall  in  a  stable  straight  line,  in  a  direction 
directly  opposite  the  pinger  location,  and  at  an  angle  with 
the  vertical  equal  or  close  to  the  package  initial  tilt 
angle. 

Initial  package  inclination  did  not  result  in  package  tumbling. 
Packages  would  simply  straighten  their  flight  back  to  vertical 
fall  and  then  resume  the  flight  pattern  specific  to  their  parti¬ 
cular  configurations.  This  seems  to  confirm  the  static  stability 
of  all  configurations  tested,  except  of  course  for  those  with 


a  nose  cone. 


In  summary,  with  the  exception  of  those  configurations  with  a  nose 
cone,  all  configurations  seem  to  be  reasonably  statically  stable.  Serious 
tumbling  or  kiting  could  not  be  detected. 

Both  the  baseline  and  the  baseline  with  a  pinger  are  dynamically  un¬ 
stable,  i.e.  they  wobble  around  their  line  of  flight.  These  two  config¬ 
urations  are  those  most  often  encountered  in  practice.  The  oscillations 
can  be  reduced  and/or  suppressed  by  either  lowering  the  c.g.  (add  weight 
close  to  the  bottom)  or  by  adding  stabilizing  fins.  The  latter  solution 
however,  will  force  the  nonsymmetrical  package  to  fall  in  a  straight  line 
and  glide  off  the  vertical.  On  the  other  hand  adding  weight  increases 
the  terminal  velocity,  a  desirable  result.  Based  on  these  results,  the  fol¬ 
lowing  recommendations  to  improve  package  configuration  and  performance 
can  be  made: 

Avoid  horizontally  mounted  instruments. 

Avoid  or  minimize  package  asymmetry. 

Judiciously  add  weight  at  the  bottom  of  the  package. 

Reduce  package  drag  by  means  other  than  a  nose  cone. 

4.0  SEA  TESTS 

4.1  Objectives 

In  order  to  confirm  the  motions  and  tension  values  predicted 
through  analyses  and  scale  model  tests  a  standard  CTD  package  was  modified 
and  taken  to  sea.  Measurements  of  line  tension  and  package  tilt  during 
routine  lowerings  were  desired  at  various  payout  speeds  and  under  different 
sea  state  conditions  to  observe  the  performance  of  the  package. 

4.2  Instrumentation 

A  standard  CTD  was  modified  to  provide  tilt  in  two  orthogonal 
axis  and  line  tension  immediately  above  the  package.  The  tilt 
meters  measure  the  angle  the  CTD  package  makes  with  the  vertical  in  order 


to  determine  whether  the  instnnnent  is  timbling  and/or  spinning.  The  meters 
were  positioned  at  right  angles  to  each  other  in  a  plane  perpendicular  to 
the  long  axis  of  the  package  nine  inches  below  the  top  end  cap  of  the  CTD. 

The  meters  are  Humphrey,  Inc.'s  critically  damped  pendulum  with  a  range  of 
±45° .  The  pendulum  drives  a  wiper  across  a  2000  ohm  potentiometer.  As 
shown  in  Figure  8a  the  tilt  sensor  outputs  a  voltage  between  0  and  5  volts. 
This  voltage  is  sensed  by  a  differential  amplifier  which  buffers  the  input 
to  a  12-bit  analog-to-digital  converter  resulting  in  a  resolution  of  .022®. 
The  combined  error  of  both  the  differential  amplifier  and  the  converter 
is  less  than  .05°  over  the  temperature  range  to  which  the  instrument  was 
exposed.  The  total  accuracy  is  better  than  .1°  for  either  tilt  sensor. 

Line  tension  was  measured  using  a  BLH  electronics  350  ohm  load  cell 
modified  to  permit  the  use  of  the  cell  in  an  ocean  environment  to  a  depth 
of  7000  meters  without  error  due  to  ambient  pressure  effects.  The  load 
cell  has  a  range  of  0  to  5800  pounds  and  infinite  resolution.  The  maximum 
error  from  all  causes  is  .2%  of  rated  output  over  the  temperature  range 
the  instrument  was  exposed  to.  As  shown  in  Figure  8b  the  output  from  the 
load  cell  is  amplified  667  times  to  give  a  0  to  10  V  output  for  a  load  of 
0  to  4000  lbs.  This  output  is  buffered  by  a  differential  amplifier  into 
a  12-bit  analog-to-digital  converter  resulting  in  a  resolution  of  one 
pound  and  a  total  system  accuracy  of  within  ±10  pounds. 

The  depth  velocity  and  acceleration  of  the  CTD  package  is  deter¬ 
mined  from  differentiation  of  the  normal  pressure  sensor  output  of  the  CTD. 

It  is  a  350  ohm  bonded  transducer  manufactured  by  Paine  Instruments.  Its  out¬ 
put  is  amplified  a  hundred  times  and  then  converted  into  a  16-bit  binary 
number.  The  CTD  package  used  has  a  pressure  range  of  0  to  6500  dbar  and  a 
resolution  of  0.1  dbar.  The  accuracy  is  ±6,5  dbar. 
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All  of  the  sensor  outputs  are  digitized  twenty-four  times  a  second 
giving  a  temporal  resolution  of  41  msec.  The  data  is  organized  in  a 
serial  data  stream  for  transmission  up  the  E/M  cable  to  a  shipboard 
pdeck  unit.  The  deck  unit  decodes  the  serial  data  for  display  and  also 
-supplies  binary  coded  audio  signals  which  are  recorded  by  a  standard  audio 

iv 

tape  recorder.  The  data  can  then  be  analyzed  at  a  later  time  by  replay- 

— . 

'•ring  the  audio  tape. 

I  .  4.3  Procedure 

The  specially  instrumented  CTD  package  was  deployed  at  sea 
*  from  the  R/V  OCEANUS  on  two  separate  cruises.  Both  lowerings  were  done 
with  the  cable  fed  through  a  sheave  attached  to  the  A-frame  on  the  fantail. 
The  normal  method  of  lowering  is  from  the  boom  on  the  starboard  side  but 
because  of  operational  limitations  on  these  cruises  the  package  was 
lowered  in  this  manner.  No  other  changes  were  made  in  the  usual  operating 
procedures.  There  was  a  major  difference  in  sea  conditions  between  the 
two  cruises.  The  first  cruise  was  during  calm  to  moderate  weather  with 
small  gentle  ship  motions.  The  second  lowering  was  done  about  six  hours 
after  an  intense  winter  storm  with  large  sharp  changes  in  pitch  and  roll 
angles.  Both  cruises  were  150  miles  south  of  Martha *s  Vineyard  in  the 
late  fall/winter  season. 

12/14/81  Cruise 

The  weather  as  reported  by  the  OCEANUS  bridge  was  calm  seas,  wind 
southwest  at  10  knots,  ship  pitch  and  roll  less  than  3®,  increasing 
slightly  after  the  CTD  reached  1700  meters  depth.  The  CTD  package  was 
attached  to  a  standard  .303  electromechanical  cable  wound  on  the  ship’s 
Markey  winch.  The  package  was  lowered  into  the  sea  and  halted  at  hull 
depth.  It  was  then  lowered  at  25  meters  per  minute  to  60  meters  where 
its  speed  was  increased  to  40  meters  per  minute.  At  228  meters  depth 


itiie  lowering  speed  was  increased  to  60  meters  per  minute  and  remained 

H '  I 

constant  until  the  bottom  of  the  cast  at  around  2000  meters  depth.  In 
order  to  simulate  a  normal  CTD  cast  rosette  bottles  were  fired  at  twenty 
Idifferent  depths  during  its  retrieval.  At  each  bottle  firing  the  CTD 

I 

Ipackage  was  stopped  for  15  to  20  seconds.  Throughout  the  cast  the  data 
Sfrom  the  CTD,  tension  cell  and  tilt  meters  was  monitored  onboard  the  ship 

Sii.-.  . 

i- 

iand  recorded  on  an  audio  tape  recorder. 

X 

?  1/6/82  Cruise 

V' 

The  sea  conditions  during  this  test  were  winds  at  10  knots  from  the 
west,  waves  at  6  feet  and  swells  also  from  the  west.  The  CTD  package  was 
lowered  into  the  water  with  a  speed  of  twenty  meters  per  minute  to  172. 
meters.  The  data  from  the  tension  cell  stopped  at  78  meters  due  to  storm 
damage  to  the  electrical  harness.  The  pressure  and  tilt  meters  continued 
to  operate.  The  CTD  package  was  retrieved  with  13  stops  for  rosette 
firings  of  the  Niskin  bottles.  The  data  was  again  recorded  on  an  audio 
tape  recorder. 

4.4  Data  Processing 

The  audio  tapes  from  the  two  lowerings  were  replayed  at 
WHOI.  The  data  was  monitored  by  a  CTD  deck  unit  with  a  RS-232  serial 
data  link  with  a  Hewlett  Packard  2113  computer.  All  of  the  data  was  trans 
cribed  onto  a  nine-track  digital  magnetic  tape  and  once-a-minute  listings 
were  made  of  each  of  the  CTD  data  channels. 

Before  and  after  each  of  the  lowerings,  the  tilt  sensors,  tension 
cell  and  pressure  transducer  were  calibrated.  These  calibrations  were 
then  used  to  determine  connection  factors  for  the  raw  data  from  the  casts. 
The  data  processing  software  allows  a  linear  correction  to  each  of  the 
data  channels.  This  corrected  data  can  then  be  displayed  as  listings 
and/or  plots  with  the  data  converted  into  engineering  units.  The  computer 
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,o<jes  can  display  every  data  point  or  every  nth  point  and  also  will 
leverage  any  number  of  observations.  Any  of  the  data  channels  can  be 
iplotted  vs  the  others  and/or  time. 

Initial  data  processing  emphasized  low-frequency  events.  Both 
orthogonal  tilt  angles  and  the  tension  were  plotted  as  a  function  of  time 

land  of  depth  once  every  second.  A  listing  of  these  data  points  was  also 

. 

liHinade.  From  these  plots  areas  of  special  interest  for  high-frequency 

ip 

prevents  were  identified.  These  areas  were  plotted  on  an  expanded  time  scale 

p:!;#-once  every  41  milliseconds.  Plots  of  pressure  vs  time  were  done  of  both  low- 

1. 

erings  for  every  data  point  and  used  to  determine  velocities  and  accelerations. 


k  4.5  Sea  Tests  Results 

An  analysis  of  the  records  of  instrument  package  depth, 


cable  tension,  and  package  inclination  was  made  to  determine  the  package 
and  cable  behavior  as  they  were  lowered  and  retrieved  through  the  water 
column. 

Events  of  particular  interest  that  this  analysis  could  help  detect 
and  quantify  included  the  following: 

o  Instantaneous  rate  of  cable  fall 

o  Instantaneous  rate  of  CTD  package  fall 

o  Cable  relaxation 

o  Cable  snap  loads 

o  Instrument  timibling  and/or  spinning. 

Entire  records  were  first  examined.  This  first  look  was  used  to  detect 
long  period  trends,  to  measure  slow  varying  parameters  -  winch  payout 
rates  for  example  -  to  determine  maxima  and  minima  of  tension  and  their  rate 
of  occurrence,  and  to  identify  record  sequences  of  special  interest.  The 
time  scale  of  these  sequences  was  later  expanded  and  the  data  studied  in 
more  detail.  The  results  of  this  analysis  are  hereafter  summarized. 

Cable  Speed  of  Fall.  Pressure  records  can  be  used  to  hindcast  the 
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^ed  of  cable  fall.  By  differentiating  the  pressure  record  one  obtains 
^record  of  instrument  rate  of  travel  in  the  vertical  direction.  On  a 
pwn  roll  the  cable  can  fall  faster  than  the  instrument  and  therefore 
down  speed  measurements  cannot  always  be  used  to  infer  cable  speed.  On  the 
other  hand  when  the  instrument  is  pulled  upwards  then  cable  and  instrument 
travel  together.  If  at  that  moment  the  winch  is  paying  the  cable  out  at 
a  given  rate  then  obviously  the  sheave  must  be  going  up  at  a  speed  equal 
to  the  sum  of  the  payout  rate  and  the  rate  of  instnnnent  climb.  Now  old 
salts  would  agree  that  in  general  a  hard  roll  is  followed  by  an  equally 
hard  down  roll.  Thus  it  seems  fair  to  "expect"  (in  a  statistical  sense) 
the  sheave  to  fall  down  as  fast  as  it  came  up.  The  cable  down  speed 
would  then  be  the  down  speed  of  the  sheave  augmented  of  course  by  the 
payout  rate,  or  the  cable's  own  terminal  velocity  whichever  is  the  smaller. 

Applying  this  reasoning  to  the  pressure  data  shown  in  Figure  9,  the 
upward  speed  of  instrument  at  time  t  =  8  seconds  is  found  to  be  5.25  ft/sec. 
The  winch  paying  out  at  a  rate  of  2,4  ft/sec,  the  sheave  must  be  climbing 
at  7.65  ft/sec.  On  the  next  down  roll  the  cable  will  fall  at  7.65+2.4  = 

10  ft/sec.  The  pressure  record  shows  the  instrument  to  fall  at  that  time 
(t  -  10  sec)  at  9.2  ft/sec.  The  cable  falling  faster  goes  slack,  as  evi¬ 
denced  by  the  tension  record.  This  example  clearly  shows  that  on  occasions 
the  cable  can  and  indeed  does  fall  faster  than  the  instrument. 

The  maximinn  upward  sheave  speed  observed  in  all  records  was  found  to 
be  11,58  ft/sec.  As  further  explained  in  Appendix  1,  when  the  shaave  falls 
at  that  rate,  its  speed  then  exceeds  the  instrument  terminal  velocity 
(8.2  ft/sec).  It  could  also  be  equal  to  or  larger  than  the  estimated  rate 
of  cable  free  fall  (12.2  ft/sec  il.2  ft/sec).  In  these  conditions,  and 
with  the  winch  secured,  the  cable  could  certainly  be  slack  at  its  lower 
end  and  also  most  probably  at  its  upper  end.  Paying  cable  out  would 
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FIGURE  9.  Records  of  tension  and  pressure  measurements 
(R/V  OCEANUS  Cruise  112) 


just  increase  the  chances  for  cable  relaxation  at  either  or  both  ends. 

Tension  Measurements .  Tension  measurements,  particularly  those  of 
the  second  lowering,  revealed  large  deviations  from  the  instrument  immersed 
weight  (306  lbs) .  Low  levels  of  tension  -  including  zero  -  were  often  en¬ 
countered  followed  by  values  as  high  as  three  times  the  instrument  weight. 

A  typical  dynamic  section  of  the  record  was  selected  for  detailed  analysis. 
vFigure  9,  based  on  24  data  points  per  second,  shows  the  synchronous  measure- 
ij.jnents  of  cable  tension  and  instrument  depth  as  a  function  of  time,  starting 
from  an  arbitrary  origin.  This  extremely  interesting  record  illustrates 
cable  relaxation  followed  by  impact  loading. 

At  time  t  =  4  seconds  the  instrument  starts  to  fall.  The  cable  over¬ 
rides  the  instrument  1.5  seconds  later  and  the  tension  goes  to  zero 
(t  =  5.5  sec).  The  cable  remains  slack  for  a  second.  At  t  =  6.5  sec  the 
instrument  starts  decelerating  and  is  brought  to  zero  speed  0.25  sec  later. 
The  tension  jumps  from  zero  to  706  lbs  over  this  short  time  interval.  The 
snap  load  is  followed  by  typical  shock  waves  in  the  cable.  As  the  instru¬ 
ment  travels  upwards  the  pull  of  the  cable  becomes  less  and  less.  Event¬ 
ually  the  instrument  stops  climbing  (t  =  9  sec),  starts  to  fall  again  ac¬ 
quiring  terminal  velocity  at  t  =  10  secs  at  which  time  the  cable  overrides 
the  instrument  again.  A  second  relaxation  takes  place,  followed  by  a  second 
snap  load  (t  =  11  sec) . 

Peak  tensions  reached  in  this  particular  series  of  events  -  706  lbs  or 
so  -  may  seem  deceivingly  small  when  compared  to  the  7300  lbs  of  cable 
strength.  On  the  other  hand  more  severe  weather  undoubtedly  would  have 
produced  much  larger  peaks.  Of  more  significance  perhaps  is  that  relaxation 
-  which  can  permit  kinks  to  form  -  happened  repeatedly.  Repeated  snap  loads 
would  them  be  applied  to  a  damaged  cable. 
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An  attempt  was  made  at  reproducing  the  cable  tension  history  using 
the  mathematical  model  described  in  Appendix  2,  As  can  be  seen  from 
Figure  9,  predictions  of  cable  relaxation  followed  by  snap  loads  agree 
rather  well  with  the  measured  values,  despite  model  simplicity. 

Tilt  Measurements.  The  records  indicated  that  the  package  was  slightly 
tilted  (5  to  6  degrees)  and  was  turning  slowly  while  being  lowered  and 
raised.  The  slow  spin  was  probably  caused  by  cable  unlaying.  Violent 
tumbling  could  not  be  observed. 

5.0  RECOMMENDATIONS 

Causes  of  lowering  cable  damage  were  identified  and  measured 
in  situ.  These  include  high  stress  levels  at  the  head  sheave,  cable 
relaxation  and  snap  loads.  The  flight  pattern  of  typical  CTD  packages  has 
been  observed  with  the  help  of  scale  models.  Math  models  were  used  to 
quantify  tension  levels,  calculate  terminal  velocities  and  investigate 
instrument  package  stability.  Based  on  this  information  the  following 
recommendations  to  reduce  or  suppress  the  cause  of  cable  damage  and  improve 
the  package  flight  pattern  are  made. 

5.1  Modus  Operand! 

As  a  first  measure,  limits  should  be  set  and  observed 
as  to  maximum  depth  of  casts  and  allowable  winch  speeds.  To  help  plan 
safe  lowerings,  predictions  of  tension  at  the  head  sheave  should  be 
readily  available.  If,  for  example,  one  had  graphs  of  peak  tension  versus 
cable  length  for  different  hauling  speeds  and  sea  states  as  shown  in 
Figure  10,.  then  the  allowable  cable  length  could  be  found  from  the  inter¬ 
section  of  the  pertinent  tension  curve  with  the  safe  load  (or  the  elastic 
limit)  of  a  particular  cable. 


CTD  lowering  operations  should  certainly  try  to  avoid  slack  cable 
conditions.  Calculations  of  winch  payout  rates  which  would  prevent  the 
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E/M  CABLE  CHARACTERtSTICS: 

IMMERSED  WEIGHT  =  .145lb/ft  (.22kg/m) 
DIAMETER  =  .303 in.  (.77cm) 

BREAKING  STRENGTH  =  7400lbs  (3360kg) 
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TYPE  :  CTD  PACKAGE 

IMMERSED  WEIGHT  =  370 lbs  (168 kgs) 


HEAVE  AMPLITUDE  =  4ft  (1.2m) 
ROLL  AMPLITUDE  =  10  DEGREES 
WAVE  PERIOD  =  6  SECONDS 


VESSEL  :  R/V  OCEANUS 


FIGURE  10.  Curves  of  peak  tension  versus  cable  length 
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cable  to  become  slack  as  a  function  of  known  (preferably  measured) 
package  terminal  velocity  and  sea  state  should  be  made  and  presented  in 
a  friendly  tabular  or  monograph  form.  Limits  on  payout  rates  should  then 
be  set  accordingly.  Examples  are  shown  in  Figure  11  of  limits  on  payout 
rates  for  two  different  packages  having  terminal  velocities  of  8  and  12 
ft/sec  respectively.  It  is  apparent  that  every  effort  should  be  made  to 
use  a  package  which  has  as  high  a  terminal  velocity  as  the  weight  limits 
will  allow. 

5.2  Improving  CTD  Packages 

Improving  existing  CTD  packages  can  be  achieved  by  de¬ 
creasing  the  projected  frontal  area  and  increasing  the  weight  by  efficient 
ballasting.  There  are  obviously  limitations  depending  on  the  size  and 
shape  of  the  original  package.  If  a  large  and  heavy  package  exists,  it 
may  be  difficult  to  increase  its  weight  efficiently  and  still  end  up  with 
one  that  can  be  practically  managed  at  sea. 

A  package  terminal  velocity  larger  than  the  terminal  velocity  of  the 
cable  would  certainly  reduce  the  chances  of  cable  slackness  and  kinking  at 
the  CTD  endo  In  the  cases  studied  here,  that  cable  velocity  is  just  over 
12  ft/sec.  To  achieve  a  velocity  equal  to  or  greater  than  this  our  standard 
CTD  package  was  redesigned.  The  mounting  frame  was  first  removed  to  reduce 
drag.  To  compensate  for  the  reduction  in  weight,  ballast  in  the  form  of 
a  lead  collar  was  attached  around  the  lower  portion  of  the  CTD  instrument 
case.  This  further  reduced  drag  by  eliminating  undesirable  end  cap  effects 
caused  by  separate  cylindrical  instnmient  cases.  The  resulting  smaller 
and  faster  package  also  accelerated  more  quickly  as  the  plots  in  Figure  6 
show.  This  efficient  configuration  can  thus  more  readily  react  to  dynamics 
caused  by  ship  motions. 
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FIGURE  11.  Curves  of  winch  payout  versus  heave  and  roll 
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I  Static  stability  is  important  and  can  be  best  maintained  by  keeping 
nstruments  fairly  long  with  heavy  items  low  (ballast)  and  light j  high 
rag  items  (bottles)  high. 

5.3  Improving  Lowering  Equipment 

Setting  limits  on  deployment  depths  to  avoid  high  stresses 
lift  the  sheave  and  on  payout  rates  to  prevent  cable  slack  conditions  should 
considered  but  temporary  measures.  Certainly  the  desirable  goal  remains 
Ithe  ability  to  perform  deeper  and  faster  casts  in  more  severe  weather 
•conditions. 

Ship  motion  produces  high  cable  stresses,  kinks,  snap  loads  and  often 
introduces  deleterious  effects  on  the  quality  of  the  data  collected. 
Automatic  motion  compensation  would  help  reduce  or  suppress  these  undesira¬ 
ble  causes  of  system  failure.  Shock  loads  can  be  substantially  reduced 
^  with  the  help  of  multiple  sheave  hydraulic  or  pneumatic  accumulators.  An 
active  boom  to  suspend  the  head  sheave  can  also  be  used.  Another  approach 
I  is  to  modulate  the  speed  of  the  winch  in  response  to  measured  sheave  dis¬ 
placement  rates.  Accumulators  and  active  booms  are  complex  systems,  they 
i;  require  ample  deck  space,  and  are  stroke  limited.  For  open  ocean  applica- 
tions,  the  servo-controlled  winch  appears  to  be  a  more  practical  and  attrac¬ 
tive  solution.  A  short  discussion  follows  of  the  benefits  introduced 
by  winch  speed  regulation. 

When  hauling  the  cable  and  the  instrument  back  on  board  ship  in  a 
rough  sea  way,  the  tension  at  the  sheave  is  the  instantaneous  sum  of  the 
steady  state  load  due  to  cable  and  instrument  immersed  weight  and  hauling 
drag  and  of  the  time  varying  load  due  to  inertia  and  drag  forces  imparted 
on  the  cable  and  the  instrument  by  the  violent  motion  of  the  head  sheave. 
These  wave  induced  forces  cause  repeated  high  tension  peaks.  In  this  in¬ 


stance  the  instantaneous  speed  of  ascent  of  the  cable  and  instrument  is 


the.  vector  sum  of  the  winch  constant  hauling  speed  and  of  the  time  depen¬ 
dent  vertical  component  of  the  sheave  speed.  If  the  winch  hauling  speed 
could  be  made  time  dependent  and  regulated  so  as  to  produce  a  constant 
vector  sum  then  the  actual  speed  of  cable  and  instrument  ascent  would  be 
constant.  Inertia  effects  resulting  from  speed  changes  would  be  suppress¬ 
ed  as  well  as  the  substantial  increase  in  drag  forces  due  to  upward  sheave 
tnotion.  The  tension  would  be  nearly  constant,  decreasing  slowly  as  the 
cable  is  reeled  in.  High  tension  peaks  due  to  dynamic  loading  would  no 
longer  be  present.  Deeper  casts  could  be  made  in  more  severe  weather 
I?  conditions . 

When  lowering  an  instrument  package  in  a  rough  sea  way  cable  slack¬ 
ness  followed  by  snap  loads  is  likely  to  occur.  Here  again  the  actual 
cable  lowering  speed  is  the  vector  sum  of  the  winch  payout  rate  and  of 
the  vertical  component  of  the  head  sheave  speed.  If  the  winch  payout 
g  rate  could  be  regulated  to  produce  a  constant  vector  sum  smaller  than  the 
instrument  terminal  velocity  then  slack  conditions,  kinks,  and  snap  loads 
would  be  entirely  eliminated.  Instruments  could  be  lowered  at  much  faster 
rates  and  in  more  severe  weather,  thus  conserving  precious  ship  time.  A 
^  scheme  to  regulate  winch  speed  is  depicted  in  Figure  12.  (Personal 
communication,  J.  G.  Desureault,  BIO,  Canada,  1982). 

To  conclude,  automatic  motion  compensation  would  greatly  reduce  the 

mik- 

causes  of  failure  identified  in  this  study.  By  enabling  casts  to  be 
made  in  larger  depths  and  more  severe  weather  it  would  put  costly  ship 
time  to  better  use.  By  providing  smooth  lowering  and  vertical  station 
gp;  keeping  capability  it  would  enhance  data  quality  and  expand  the  scope  of 
scientific  experiments  using  profiling  instrumentation  techniques. 

The  need  to  improve  and  upgrade  present  oceanographic  winches 
has  been  fully  recognized  and  motion  compensation  systems  best  suited 
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FIGURE  12.  Servo  controlled  winch  schematic 


^PPEM)IX  1»  Maximim  Sheave  Speed  and  Estimated  Terminal  Velocity  of 


CTD  Cable 
Maximum  Sheave  Speed 


Figure  13 

The  section  of  the  GTD  depth  record  with  the  largest  recorded 
rate  of  instrument  ascent  is  shown  in  Figure  13  above.  From 
this  record,  the  rate  is  found  to  be  2.8  m/sec  or  9.18  ft/ sec. 
Adding  2.4  ft/sec  for  the  payout  rate  at  the  time  yields  a 
measured  sheave  upward  speed  of  11.58  ft/ sec. 
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Cable  Terminal  Velocity 

The  highest  possible  cable  free  fall  velocity  is  achieved  when  the 
cable  has  a  minimum  drag-to-weight  ratio.  This  obviously  is  ob¬ 
tained  when  the  cable  falls  vertically. 

Let  us  consider  a  length  of  cable  L  falling  vertically.  At 
terminal  velocity  Vrp,  the  gravity  pull  equals  the  cable  drag. 

This  condition  can  be  expressed  by: 

where  =  -r 

=  Cable  immersed  unit  weight  (Ib/dt) 

^  =  Sea  water  mass  density  =  2  slugs/cu.ft 

Q  =  Cable  diameter  (ft) 

^ =  Gable  longitudinal  drag  coefficient. 

From  this  equality  the  terminal  velocity  is: 


Using  values  of 

D  =  0.303  in  (CTD  cable  specs) 

\A/^^  =  0.116  Ib/dt  (CTD  cable  specs) 

=  .012  (high)  and  .008  (low) 

yie Ids 

High  terminal  velocity  =  13.52  ft/sec  =  (4.03  m/sec) 

Low  terminal  velocity  =  11.03  ft/sec  =  (3.36  m/sec) 

Average  terminal  velocity  =  12.27  ft/sec  =  (3.70  m/sec) 

Thus  the  estimated  speed  of  maximum  cable  free  fall  is 
12.27  ft/sec  (±  1.25  ft/sec) 

or  3.74  m.sec  (±  0.34  m/sec) 

or  224  m/min  (±  20  m/min) 


for  the  need  are  being  studied. 

A  last  remark  should  be  made  to  emphasize  once  more  the  importance 
of  good  cable  handling  practices.  A  wealth  of  information  on  this 
critical  subject  can  be  found  in  the  recently  published  Handbook  on 
Oceanographic  Winch,  Wire  and  Cable  Technology  (Reference  4), 
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I;  appendix  2.  Mathematical  Model  of  CTD  Package  Dynamics 

•.’A  , 

In  the  course  of  the  CTD  study  reported  herein  an  attempt  was  made 
at  developing  and  exercising  a  simple  computer  model  which  could  predict 
the  occurrrence  of  cable  relaxation  followed  by  snap  loads. 

Following  the  approach  described  in  Reference  5,  the  single  degree 
of  freedom  system  shown  in  Figure  14  is  used  to  represent  the  package 
and  lowering  cable.  In  this  model  the  following  assumptions  are  made: 

o  The  motion  of  the  package  is  entirely  vertical. 

o  The  mass  of  the  cable  is  assumed  to  be  a  small  fraction  of  the 
equipment  mass.  This  would  be  the  case  for  rather  short 
lengths  of  cable  (hundreds  of  meters  instead  of  thousands) , 
or  if  the  cable  is  light  (Kevlar  line  for  example),  or  if  the 
payload  entrains  a  lot  of  water. 

o  The  cable  is  treated  as  an  elastic  spring  which  changes  its 
length  as  the  winch  pays  the  cable  out  or  reels  it  in. 

In  the  elastic  range  of  cable  elongation  the  spring  constant  k 
is  given  by 


K  = 


where  £  is  the  cable  modulus  of  elasticity  (psi) 
A  is  the  cable  metallic  area  (sq/in) 

U  is  the  cable  length  (ft) • 

K,  will  then  be  expressed  in  Ibs.ft, 
Equipment  drag  is  nonlinear  and  of  the  form 


o  Equipment  drag  is  assumed  much  larger  than  the  cable  drag. 

The  cable  upper  end  (Xj(*f))  is  forced  to  follow  a  prescribed  dis¬ 
placement  such  as  a  sinusoidal  or  a  known  time  series.  The  ensuing  dis¬ 
placement  of  the  CTD  package  ( (i)  )  is  found  by  numerical  integration 
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FIGURE  14.  Mathematical  model  of  CTD/cable  system 
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:!•  the  instrument  equation  of  motion,  namely: 

K  ( -  Xg)  ■“  Ap  I  1  ■=  lOT)^  X , 


(2.1) 


where 


I ..  ^  is  the  sea  water  mass  density  (slug/cu.ft) 

,■ _  is  the  drag  coefficient  of  the  package 

T  Ap  is  the  cross  section  area  of  the  package  normal 

to  the  direction  of  travel  (sq.ft) 
is  the  speed  of  package  travel  (ft/sec) 
is  the  package  acceleration 
is  the  package  virtual  mass. 

The  instantaneous  cable  tension  at  the  instrument  is  in  turn  found  from 

(2.2) 

V/p  being  the  package  immersed  weight 
When  the  tension  goes  to  zero,  the  motion  of  the  instrument  is 

governed  by  the  "free  flight"  equation: 


-  Wp-  j 


(2.3) 


When  the  tension  becomes  positive,  Equation  (2.1)  prevails  again.  Speed  and 
displacements  used  when  switching  from  one  equation  to  the  other  are  of 
course  those  computed  at  the  time  immediately  preceding  the  switch  over. 

A  FORTRAN  computer  program  (GTDDOWN)  was  written  (see  listing  at 
end  of  this  Appendix)  to  numerically  integrate  the  equation  of  motion. 
Typical  program  input  includes  type  of  surface  displacement  prescribed, 
payout  rates,  cable  characteristics  and  instrument  package  characteristics. 
Typical  program  output  includes  a  summary  of  the  input  data,  and  a  listing 
of  head  sheave  displacement  and  speed,  of  package  depth  and  speed,  and  of 
cable  tension  as  a  function  of  time. 


Page  3  of  A. 2 


-This  computer  program  was  used  to  investigate  the  behavior  of 

If'/- 

icltl  CTD  packages  when  deployed  from  a  rolling  and  heaving  ship  of 
X  geometry.  Using  a  suitable  time  series  of  sheave  displacement, 

-ues  of  cable  tension  and  GTD  depth  were  computed  and  compared  against 
ual  measurements.  When  using  a  spring  constant  K  as  previously 
ned,  the  tensions  predicted  were  found  to  be  larger  than  those  recorded, 
ver  reducing  the  stiffness  of  the  theoretical  constant  did  produce 
reasonable  predictions  shown  in  Figure  9. 
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TYPE  CTDD0WN»F0RJ28 

CTDDOWN  PROGRAM. THIS  PROGRAM  COMPUTES  THE  INSTANTANEOUS  DEPTH  AND 
SPEED  OF  A  CTD  INSTRUMENT  PACKAGE  BEING  LOWERED  FORM  A  ROLLING  SHIP. 
THE  TENSION  AT  THE  CABLE  CTD  END  IS  ALSO  COMPUTED. 

PROGRAM  WRITTEN  BY  H.O.BERTEAUX, JUNE  1  1983 

PROGRAM  REOUIRED  INPUT  INCLUDES? 

SPECIFICATION  OF  SURFACE  EXCITATION 

CHARACTERISTICS  OF  LOWERING  CABLE, AND  LOWERED  PACKAGE 
INITIAL  CONDITIONS 

NON  EXECUTABLE  STATEMENTS 

i 

CHARACTERJKSO  FILENAME ,  CAbLeNAME ,  INSTRUMENT ,  SHIPNAME , DATE 
DIMENSION  SHEAUEdOOO)  ,TIMDAT(1000> 


INPUT  STATEMENTS 
IDENTIFIES  SHIP  AND  DATE 
WRITE ( 6, ))() 'NAME  OF  SHIP?' 

READ<5rl000)  SHIPNAME 
FORMAT (A) 

WRITE<6,)|c)  'DATE?' 

READ (5 ,1000)  DATE 

SPECIFIES  SURFACE  EXCITATION 

WRITE<6,)K) 'IS  SHEAVE  DISPLACEMENT  SPECIFIED  IN  DATA  FILE'?' 
WRITE<6,)(C) 'IF  YES  ENTER  l^^.IF  NO  ENTER  0' 

READ<5,)»:)  ANSWER 
IF  (ANSWER^EQd)  THEN 
GO  TO  200 
ELSE 

GO  TO  100 
END  IF 

WRITE ( 6, )|c) 'SHEAVE  DISPLACEMENT  IS  ASSUMED  TO  BE  PERIODIC' 
WRITE<6,)K) 'ENTER  HEAVE  AMPLITUDE < FT )  ' 

READCSrJfc)  HEAVA 

WRITE < 6, «) 'ENTER  ROLL  AMPLITUDE < DEGREES)  AND  PERIOD ( SECS ) ' 
READ<5,)»c)  ROLLA,ROLLT 

WRITE(6,)K) 'ENTER  HORIZONTAL  AND  VERTICAL  DISTANCE  TO  SHFAVFCFT)' 
READ(5,)Jc)  XSHEAVE,YSHEAVE 

WRITE(6,)(c) 'ENTER  COMPUTATION  TIME  INCREMENT ( SECS ) ' 

READ<5,>lc)  TINC 
GO  TO  400 

WRITE (6,*) 'NAME  OF  DATA  FILE?' 

READ (5, 1000)  FILENAME 

WRITE ( 6, ){c) 'ENTER  TIME  SPAN  AND  TIME  INTERVAL  BETWEEN  DATA  POINTS' 
READ<5,)ic)  TSPAN,TDATA 

WRITE < 6, >{c) 'ENTER  COMPUTATION  TIME  INCREMENT' 

READ<5,)(c)  TINC 

WRITE < 6, :«c) 'ENTER  HORIZONTAL  AND  VERTICAL  DISTANCE 
TO  SHEAVE (FT)' 

READ(5,)lc)  XSHEAVEfYSHEAVE 
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IP0INTS=<TSPAN/TDATA)+1 

OPEN ( UNIT=1 , FILE=FILENAHE , STATUS= ' OLD ' , READONLY ) 

DO  300  I=1,IP0INTS 
READ<1,)»:)  TIMDAT(I)  ySHEAVEd) 

SHEAOE  <  I )  =:SHEAOE  <  f)  >IC3  ♦  281 
CONTINUE 

END  OF  SURFACE  EXCITATION  INPUT 


NOW  ENTERS  INSTRUMENT  PACKAGE  CHARACTERISTICS 

WRITE ( 6, )|c) 'ENTER  CTD  CHARACTERISTICS' 

WRITE  <6r)Jc) 'ENTER  INSTRUMENT  NAME' 

READ <5r 1000)  INSTRUMENT 

WRITE<6r:|c) 'ENTER  CTD  VIRTUAL  MASS<SLUGS)  AND  CTD  IMMERSED  WEIGHT(IBS)' 
READ  VMASSfWETW 

WRITE<6,JK) 'ENTER  CTD  DRAG  CONSTANT < CD>{cA , SQ-FT ) ' 

READ<5,)|c)  DRAGC 

NOW  ENTERS  CABLE  CHARACTERISTICS 

WRITE(6,)|c) 'ENTER  CABLE  CHARACTERISTSICS' 

WRITE(6,)|c) 'ENTER  CABLE  NAME' 

READ(5,1000)  CABLENAME 

WRITE < 6,  JK) 'ENTER  CABLE  ELASTIC  CONSTSANT ( E)|CA ) ' 

READ<5,)lc)  ELASTC 


NOW  ENTERS  INITIAL  CONDITIONS 
WRITE <6r)|c) 'ENTER  INITIAL  CONDITIONS' 

WRITE <6r») 'ENTER  CTD  DEPTH<M)  AND  CTD  SPEED < M/SEC)  AT  TIME  ZERO' 
READ<5r)|c)  TZDrTZV 

WRITE(6f)Jc) 'ENTER  CABLE  TENSION(LB>  AND  PAYOUT  RATE(M/MIN)  AT  TIME  ZERO' 
READ<5,)|c)  TZTEN,TZPOUR 


END  OF  INPUT  DATA 

WRITE<6f  *) '*)|c)|C)lc>K  END  OF  INPUT  DATA 
NOW  PRINTS  A  SUMMARY  OF  INPUT  DATA 

WRITE<6rl010) 

F0RMAT<///,10X, 'SUMMARY  OF  INPUT  DATA',///) 

WRITE(6,1015)  DATE 

FORMAT  <  2X , ' DATE ' , 2X , A20 , // ) 

WRITE (6, 1020)  SHIPNAME 

FORMAT <2X, 'SHIP  NAME: ' ,2X, A20,//) 

IF  <ANSWER*EQ*1)  GO  TO  550 

WRITE<6r)K) 'SHEAVE  DISPLACEMENT  IS  ASSUMED  TO  BE  PERIODIC' 
WRITE<6,1025)  HEAVA 

F0RMAT<2X, 'AMPLITUDE  OF  HEAVE=' , IX ,F4 • 1 , IX , ' < FT ) ' ) 
WRITE<6,1030)  ROLLA 

FORMAT <2X, 'AMPLITUDE  OF  ROLL=' , IX, F4* 1 , IX, ' (DEGREES) ' ) 
WRITE<6,1035)  ROLLT 

FORMAT ( 2X , ' PER 1 0D= ' , IX , F4  *  1 , 1 X , ' <  SEC ) ' ) 

CONTINUE 

WRITE (6, 1040)  XSHEAVE 

F0RMAT<2X, 'HORIZONTAL  DISTANCE  FROM  SHIP  C>L* 

TO  SHEAVE=',1X,F6*2,1X,'(FT)') 
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WRITE<6,1045)  YSHEAOE 

F0RMAT<2X, 'VERTICAL  DISTANCE  FROM  SHIP 
TO  SHEAVE=' ,1X,F6*2,1X, ' (FT) ',///) 

IF (ANSWER  tEQ*  0)  GO  TO  600 
WRITE(6,1050)  FILENAME 

FORMAT (2X' SHEAVE  DISPLACEMENT  AS  SPECIFIED  IN  FILE' , IX, A50,///) 
NOW  PRINTS  CABLE  DATA 
WRITE (6,1052) 

F0RMAT(2X, 'CABLE  CHARACTERISTICS' ,/) 

WRITE(6,1055)  CABLENAME 
F0RMAT(2X, 'CABLE  NAME t ' ,2X, A20) 

WRITE(6,1060)  ELASTC 

FORMAT (2X, 'CABLE  ELASTIC  CONSTANT=' ,1X,F9»0,///) 


NOW  PRINTS  CTD  CHARACTERISTICS 
WRITE(6,1062) 

FORMAT ( 2X , ' CTD  CHARACTERI ST ICS',/) 

WRITE  (6,1065)  INSTRUMENT 

FORMAT ( 2X , ' I NSTRUMENT  NAME  t ' , 2X , A20 ) 

WRITE(6,1070)  VMASS 

FORMAT (2X, 'CTD  VIRTUAL  MASS=' ,1X,F7*2,1X, ' (SLUGS) ' ) 
WRITE(6,1075)  WETW 

FORMAT (2X, 'CTD  WEIGHT  IN  WATER' ,1X,F7,2,1X, ' (LBS) ' ) 
WRITE(6,1080)  DRAGC 

F0RMAT(2X,'CTD  DRAG  CONSTANT^ ' , IX ,F6 ♦ 2 , IX , ' ( SQ-FT ) ' , /// ) 


NOW  PRINTS  CONDITIONS  AT  TIME  ZERO 
WRITE(6,1082) 

F0RMAT(2X, 'INITIAL  CONDITIONS',/) 

WRITE (6, 1085)  TZD 

FORMAT (2X, 'CTD  INITIAL  DEPTH=' , 1X,F8^2, IX, ' (M) ' ) 
WRITE(6,1090)  TZV 

FORMAT (2X, 'CTD  INITIAL  SPEED=^' ,  1X,F6*2,1X, '  (M/SEC)  '  ) 
WRITE(6,1095)  TZTEN 

FORMAT (2X, 'CABLE  INITIAL  TENSION^' , 1X,F8.2, IX, ' (LBS) ' ) 
WRITE (6, 1100)  TZPOUR 

F0RMAT(2X, 'CABLE  PAYOUT  RATE=' , 1X,F6*2, IX, ' (M/MIN) ' ,///) 


NOW  PRINTS  HEADINGS  OF  COMPUTATIONS  RESULTS 
WRITE(6,1105) 

F0RMAT(24X, 'COMPUTATION  RESULTS' ,//) 

WRITE(6,1110) 

FORMAT (2X, 'TIME' ,3X ,' SHEAVE ', 6X ,' CTD  DEPTH' ,4X, 

'CTD  SPEED' ,2X, 'CTD  ACCEL *', 3X ,' CABLE  TENSION') 
WRITE(6,1115) 

F0RMAT(2X, ' (SEC) ' ,2X, ' (FEET) ' , 3X ,'( FEET )', 3X ,'( METER ) ' , 
2X, ' (FT/SEC) ' ,2X, ' (FT/SEC-SQ) ' ,6X, ' (LBS) ' ,/) 


NOW  PROCEEDS  WITH  ANALYSIS  OF  CTD  BEHAVIOR 
IF  (ANSWER^EQ.l)  GO  TO  2000 

COMPUTES  CTD  RESPONSE  TO  PERIODIC  EXCITATION 

PRELIMINARY  COMPUTATIONS 

COMPUTES  INITIAL  SPRING  CONSTANT 

CONVERTS  ROLL  AMPLITUDE  FROM  DEGREES  TO  RADIANS 
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R0LLA=R0LLA)fc3»  1416/180*0 

RADIUS=SQRT(XSHEAOE>>:)»c2+YSHEAOE3|c>Jc2) 
BETHA= ATAN ( YSHEAOE/XSHEAOE ) 
SURFX=RADIUS)icSIN<BETHA)-HEAOA 
DEPTH=~TZD>IC3*281 
SPRING=:ELASTC/  ( SURFX-DEPTH ) 


COMPUTES  UNSTRECHED  LENGTH  OF  CTD  CABLE 

RELAXL=SURFX~DEPTH-TZTEN/SPRING 

COMPUTES  INITIAL  CTD  SPEED 

PAY0UT=-3  *  281 5»cTZP0UR/60  *  0 

SPEED=TZV)lc3*281 

TENSION=TZTEN 

COMPUTES  INITIAL  CTD  ACCELERATION 

ACCEL=  <  TENSI ON-WETW-DRAGC))cSPEEri)fcABS  <  SPEED ) )  /OMASS 

IF <ABS( ACCEL)  ♦LE*  0*01)  ACCEL=0*0 

TIME=0*0 

CALL  RESULTS ( TIME , SURFX , DEPTH , SPEED y ACCEL y TENSION ) 


STARTS  INTEGRATION  PROCESS 
K=1 

CONTINUE 

TIME=K)|cTINC 

IF<TIME*GT*ROLLT)  GO  TO  5000 

OLDSURFX-SURFX 

OLDDEPTH=DEPTH 

SURFX=~HEAOA)»cCOS  <  2  *  0:tc3  *  14163KTIME/ROLLT )  +RADIUS5»c 
SIN<BETHA+R0LLA)KSIN<2*0}|c3*1416){cTIME/R0LLT)  ) 

CALL  INTEGRATION < ELASTC y RELAXL y PAYOUT , TIME y SPEED y 
ACCEL  y  TI NC  y  DEPTH  y  SURFX  y  OLDSURFX  y  OLDDEPTH  y  TENSION  y 
WETWyDRAGCyOMASS) 


CALL  RESULTS < TIME y SURFX y DEPTH y SPEED y ACCEL y TENSION ) 
K=K+1 

GO  TO  1500 


PROGRAM  COMPUTES  CTD  RESPONSE  TO  SPECIFIED  SHEAVE  DISPLACEMENT 
PRELIMINARY  COMPUTATIONS 

COMPUTES  INITIAL  CABLE  SPRING  CONSTANT  (SPRING) 
RADIUS=SQRT<XSHEAVEXc>[c2+YSHEAVE>fc>K2) 

BETH A=AT AN  <  YSHE A VE/XSHE A VE ) 

SURFX=SHEAVE  ( 1 )  +RADIUS:{cSIN  <  BETHA ) 

DEPTH=~TZD)ic3*281 
SPRING=ELASTC/ ( SURFX-DEPTH ) 

COMPUTES  UNSTRETCHED  LENGTH  OF  CTD  CABLE 

RELAXL=SURFX-DEPTH-TZTEN/SPRING 

COMPUTES  INITIAL  CTD  SPEED 

SPEED=TZVHC3*281 

TENSION=TZTEN 

COMPUTES  INITIAL  CTD  ACCELERATION 
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ACCEL=  ( TENSI  0N~WETW~riRAGC5tcSPEEri)|cABS  <  SPEED ) )  /VMASS 
TIME=:0*0 

CALL  RESULTS (TIME, SURFX , DEPTH , SPEED r ACCEL , TENSION ) 
COMPUTES  PAYOUTRATE 
PAY0UT=-3  ♦  281>KTZP0UR/60  ♦  0 


STARTS  INTEGRATION  PROCESS 

K=1 

M=:2 

CONTINUE 

TIME^KHcTINC 

IF  (TIME*GT*TSPAN)  GO  TO  5000 

OLDSURFX=SURFX 

OLDDEPTH=DEPTH 

CALL  DATASHEAUE < M , TIMDAT , SHEAVE , TI ME , SURFX y TINC , TDATA ) 

CALL  INTEGRATION ( ELASTC , RELAXL , PAYOUT , T IME , SPEED , ACCEL r 
TINC  f  DEPTH , SURFX , OLDSURFX , OLDDEPTH , TENSION , WETW , DRAGC , 
VMASS) 

CALL  RESULTS ( TIME , SURFX , DEPTH r SPEED , ACCEL , TENSION > 

K»K+1 

GO  TO  2100 


STOP 


END 


SUBROUTINE  D ATASHEAVE < M , TI MDAT , SHEAVE r TIME , SURFX , T I NC , 
TDATA) 

THIS  SUBROUTINE  COMPUTES  THE  SHEAVE  POSITION  (SURFX) 

FOR  EACH  TIME  INCREMENT  BY  LINEAR  INTERPOLATION  BETWEEN 
GIVEN  SHEAVE  DISPLACEMENT  DATA  POINTS 


DIMENSION  TIMDAT< 1000) ,SHEAVE( 1000) 

IF(TIME*GT»TIMDAT(M))  GO  TO  3100 
SURFX=-SURFX+  ( SHEAVE  ( M ) -SHEAVE ( M-1 )  ) )KTINC/TDATA 
GO  TO  3200 
M=M+1 

GO  TO  3000 
RETURN 


END 


SUBROUTINE  INTEGRATION<ELASTC,RELAXL,PAYOUTrTIME,SPFED, 
ACCEL , TINC , DEPTH , SURFX , OLDSURFX , OLDDEPTH , TENSION , WETW , 
DRAGCrVMASS) 
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THIS  SUBROUTINE  COMPUTES  FOR  EACH  TIME  INCREMENT  THE  CTD 
DEPTH, THE  CTD  SPEED, THE  CTD  ACCELERATION, AND  THE  CABLE 
TENSION  AT  THE  CTD  END  OF  THE  CABLE 


SPRING=ELASTC/<RELAXL~PAYOUT5|cTIME) 

SPEED=SPEED+ACCEL)|cTINC 

DEPTH=OLDDEPTH+SPEED5KTINC 

DELTALENGTH= ( SURFX-DEPTH ) - ( OLDSURFX-OLDDEPTH ) 
DELTALENGTH=DELTALENGTH+PAYOUT:»cTINC 
IF(ABS<DELTALENGTH)  ♦LE*  *001)  DELTALEN6TH=0 ♦ 0 
TENSION=TENSION+SPRING)»cDELTALENGTH 
IF(TENSI0N*LE*0*01)  TENSI0N=0*0 

ACCEL=  < TENSI0N-WETW-DRAGC3JcSPEED)|cABS <  SPEED ) ) /VMASS 
I F  <  ABS  <  ACCEL ) *  LE  *  0  *  00 1 )  ACCEL=0  *  0 


RETURN 

END 


SUBROUTINE  RESULTS < TIME , SURFX , DEPTH , SPEED , ACCEL , TENSION ) 
THIS  SUBROUTINE  PRINTS, SHEAVE  DISPLACEMENT , CTD  DEPTH, 

CTD  SPEED, CTD  ACCELERATION , AND  CABLE  TENSION  FOR  EACH 
TIME  INCREMENT* 

DEPTH1=-DEPTH 
DEPTH2=-DEPTH/3  *  281 

WRITE (6,1200)  TIME , SURFX , DEPTHl , DEPTH2 , SPEED , ACCEL , TENSI ON 
FORMAT < 1X,F6*2,2X,F6*2,2X,F8*2,2X,F7*2,2X,F7*2, 
5X,F7*2,6X,F8*2) 


RETURN 


END 


